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The tensor product representation of quantum states leads to a promising variational approach to 
study quantum phase and quantum phase transitions, especially topological ordered phases which 
are impossible to handle with conventional methods due to their long range entanglement. However, 
an important issue arises when we use tensor product states (TPS) as variational states to find the 
ground state of a Hamiltonian: can arbitrary variations in the tensors that represent ground state of 
a Hamiltonian be induced by local perturbations to the Hamiltonian? Starting from a tensor product 
state which is the exact ground state of a Hamiltonian with Z2 topological order, we show that, 
surprisingly, not all variations of the tensors correspond to the variation of the ground state caused 
by local perturbations of the Hamiltonian. Even in the absence of any symmetry requirement of the 
perturbed Hamiltonian, one necessary condition for the variations of the tensors to be physical is 
that they respect certain Z2 symmetry. We support this claim by calculating explicitly the change in 
topological entanglement entropy with different variations in the tensors. This finding will provide 
important guidance to numerical variational study of topological phase and phase transitions. It is 
also a crucial step in using TPS to study universal properties of a quantum phase and its topological 
order. 

PACS numbers: 



I. INTRODUCTION 

The central task in the study of quantum many-body 
systems is the classification of possible phases of mat- 
ter and the understanding of phase transitions between 
them. Of particular interest is the study of systems at 
zero temperature, where a whole new realm of quantum 
effects emerge, and that is what we will focus on in this 
paper. Landau's general principle for understanding con- 
tinuous phase transitions based on symmetry breaking 
and local order parameter does not apply to all phases 
and phase transitions. Topological orderj^in particular, 
is not related to any symmetry properties and topological 
phase transitions may happen between systems with the 
same^H^ or incompatible^ symmetries. It has been the 
subject of intensive research and will be the central topic 
of this paper. Aside from the lack of a qualitative under- 
standing, what makes the problem harder is the fact that 
for quantum systems whose components strongly interact 
with each other, direct numerical approach is limited due 
to strong entanglement among the particles. Generically, 
the space required for description of a quantum system 
grows exponentially with system size, hence limiting di- 
rect numerical simulations to systems usually too small 
for any practical purpose. 

Recently, insights into quantum many-body systems 
from both condensed matter physics and quantum in- 
formation science have led to the discovery of the ten- 
sor product representation of quantum states (also called 
the projected entangled pair states pH12l j which provides 
a promising variational approach to study zero tempera- 
ture quantum phase and phase transitions. Representing 



quantum many-body states with a network of tensors, 
tensor product states (TPS) are proven to be efficient for 
the study of one dimensional quantum systemJSHHI. The 
higher dimensional generalization of this approach may 
not be as efficient, yet study has shown that it repro- 
duces many known results and may reveal new features in 
systems not solvable in any conventional wajtf^HlSl The 
strength of the approach lies in the fact TPS can describe 
long range entanglement that are p resent in a large class 
of topologically ordered statesP2E3 So the variational ap- 
proach based on TPS can include both topologically or- 
dered states and symmetry breaking states and can pro- 
duce a phase diagram that contains both types of states. 
In contrast, the conventional mean-field/variational ap- 
proaches are based on states with no long range entan- 
glement, which exclude the topologically ordered states 
from the very beginning. We also note that the entan- 
glement of a simple TPS satisfies an area law^, which 
coincides with the scaling of entanglement in the ground 
state of most known non-critical systems^. 

In the variational approach based on TPS, we try to 
find a TPS which minimizes the average energy of a local 
Hamiltonian. As we change the Hamiltonian by adding 
perturbation, the tensors in the TPS are also changed in 
order to minimize the average energy for the new Hamil- 
tonian. While local physical perturbations can always be 
reflected by variations in the tensors, the other direction 
of this problem remains unclear: can an arbitrary varia- 
tion of the tensor be induced by a local perturbation of 
the Hamiltonian? 

This is a very important question if we want to dis- 
cuss phase based on states. Because phase is defined as 
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a region in Hamiltonian space, where any two points H±, 
Hi within the region can be connected by a smooth path 
without encounter singularities (i.e. phase transitions). 
So the question about tensor and phase becomes, which 
set of states in Hilbert space correspond to such a region 
in Hamiltonian space and which set of tensors in the ten- 
sor space represent these states. Starting from one point 
in the phase region, we would like to know what kind 
of variations in the tensors correspond to local perturba- 
tions to the Hamiltonian. 

We can discuss this important question in more con- 
crete setting. Assume a TPS *S?t minimizing the average 
energy of a Hamiltonian H has a property. We would like 
to ask if the property is a universal property of a phase, or 
just a special property of H . If the property is a universal 
property of a phase, then the ground state + A^t for 
the perturbed Hamiltonian H + AH still has the same 
property. If the property is a special property of H, the 
ground state for the perturbed Hamiltonian will lose this 
property. It is the collection of universal properties that 
defines a phase. So a study of universal properties is a 
study of phases. If all the variations of the tensors can 
be induced by local perturbations of the Hamiltonian, 
then we can study the stability of a property against lo- 
cal perturbation AH of the Hamiltonian by studying the 
stability of a property against variations of the tensors. 
This will give us a powerful tool to study phases using 
TPS. 

Unfortunately, it turns out that not all variations of 
the tensors can be induced by local perturbations of the 
Hamiltonian, as we show in this paper. For a generic 
TPS, which satisfies a condition called injectivitySH, ten- 
sor variations indeed correspond to Hamiltonian pertur- 
bations. However, this is not true in the general case, as 
we show in this paper with a special system with topolog- 
ical order. So it is not easy to study universal properties 
and phases using TPS. In order to use TPS to study 
phases and phases transitions, we need to find the subset 
of variations in tensors that are physical, i.e. correspond- 
ing to local perturbations of the Hamiltonian. 

For clarity, we will always refer to small changes in the 
Hamiltonian as 'perturbation' and to those in the tensors 
as 'variation'. Without any efficient method to solve for 
exact TPS representation of ground states of quantum 
many-body systems, finding the subset of the variations 
of the tensors that can be induced by local perturbation 
of Hamiltonian is in general very difficult. 

We want to, in particular, study this problem for topo- 
logically ordered phases. As TPS can give a simple de- 
scription of a large class of topological ordered states, we 
expect that it might provide a powerful tool for studying 
topological phases in general. As we know, topologically 
ordered phases are proven to be st able ag ainst any lo- 
cal perturbations of the Hamiltonia n 1 3 * 25 * 2 -^. That is, the 
topological properties, such as gr ound state degeneracy^ 
and quasi-particle statistic^ 27 * 281 , are robust under any 
local perturbation to the Hamiltonian. So in the TPS 
approach to topologically ordered phase, it is natural to 



ask: are those topological properties robust against any 
variation of tensors, that is, for any tensor which rep- 
resents a topologically ordered state, is the topological 
order robust against arbitrary variation in the tensors? 
Surprisingly, we find that this is not true. 

We focus on the Z 2 topological order represented by an 
ideal TPS in this paper and study how the topological or- 
der of the state changes as we vary certain parameters in 
the representing tensors. We characterize topological or- 
der bycalculating the topological entanglement entropy 
c^| 29 | 30 | f or s t a t e anc j observe that topological order 
(i.e. the topological entanglement entropy St p ) is sta- 
ble only against variations of the tensors that preserve 
certain Z 2 symmetry of the tensors. Since the topolog- 
ical order is robust against any local perturbations of 
Hamiltonian, this result shows that not all variations of 
the tensors correspond to local perturbations of Hamilto- 
nian. For this Z 2 model, we show that in the generic case 
Z 2 symmetry is a necessary condition for the variations 
in tensors to correspond to physical perturbations of the 
Hamiltonian. This claim is further supported by checking 
stability of the topological Renyi entropy of TPS with Z 2 
symmetry preserving variations and Z 2 symmetry break- 
ing variations of the tensors respectively. 

While calculating S tp for a general state is exponen- 
tially harcP^, we find efficient ways to do so for two sets 
of TPS near the ideal TPS with Z 2 topological order. 
For a general TPS, we calculate the topological Renyi 
entropy by mapping it to the contraction of a 2D tensor 
network, which is accomplished by using the Tensor En- 
tanglement Renormalization Algorithm^. Hence we are 
able to calculate topological entropy for regions much 
larger than was possible previously and determine the 
topological order of the state more accurately. 

Our result on the stability of topological order will help 
us in the TPS based variational approach to Z 2 topolog- 
ical phase: we should only consider the variations of the 
tensors within the subspace of tensors with Z 2 symme- 
try. The Z 2 symmetry condition and possibly other con- 
ditions will help us to understand the physical variations 
of tensors in TPS. This is crucial in using TPS to study 
quantum phases and quantum phase transitions. It may 
even lead to a classification of topological order. 

This paper is organized as follows. We start by in- 
troducing an 'ideal' lattice spin model with Z 2 topolog- 
ical order and show how the presence of topological or- 
der in the ground state wave function can be understood 
nicely with a physical mechanism called 'string-net con- 
densation'. Such a physical picture naturally gives rise 
to a simple tensor product representation of the wave 
function, to which we then add two kinds of local vari- 
ations, 'string tension' and 'end of strings'. By calcu- 
lating topological entanglement entropy numerically for 
the first case and analytically for the second case, we 
show how topological order is stable against Z 2 preserv- 
ing variations ('string tension'), but breaks down imme- 
diately when Z 2 symmetry is broken (by 'end of strings'). 
We then randomly picked 200 tensors in the neighbor- 
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hood of the ideal Z 2 TPS and calculate the topological 
Renyi entropy of the corresponding states. Tensors with 
and without Z 2 symmetry demonstrate totally different 
behavior as system size scales up. We discuss in the 
last section the implications of our findings in variational 
studies of topological phase and phase transitions. The 
details of the calculations are given in the appendix. 



II. MODELS AND RESULT 
A. Spin model with Z2 topological order 

We start from an exactly solvable model which has 
Z 2 topological ordei^^U. In this section, we give the 
system Hamiltonian, find the ground state wave func- 
tion and explain its structure and how that leads to a 
nontrivial topological order which can be detected with 
topological entanglement entropy. With these insights 
about the state we then present a simple tensor product 
representation of this wave function. 

The model is defined on a two-dimensional hexagonal 
lattice where each link is occupied by a qubit (spin-1/2). 
The Hamiltonian is a sum of commuting projection op- 
erators 



II- 



(i) 



p iep 



v j£v 



X and Z are qubit Pauli operators defined as X = 
^q^,Z=^q ^ ~\ . p stands for each hexagon pla- 

quette in the lattice and Yliep Xi is the tensor product of 
six X operators around the plaquette. v stands for each 
vertex and Yijev is the tensor product of three Z op- 
erators connected to the vertex. The ground state wave 
function has a nice interpretation using the 'string-net' 
picture where state |0) corresponds to no string on a link 
and state |1) corresponds to the presence of a string. The 
vertex term Jljeu ^3 enforces that there are even number 
of strings connected to each vertex and hence the strings 
form closed loop while the plaquette term Yl iep Xi gives 
dynamics to the closed loops. The ground state wave 
function is an equal weight superposition of all closed 
loop configurations on the lattice. 



1$. 



Ei^> 



(2) 



The normalization factor is omitted. If we refer to each 
closed loop configuration as a string-net, the appearance 
of Z 2 topological order in this system has then a natural 
interpretation as being due to the condensation of string- 
nets. We will refer to this model as the ideal Z 2 model. 

For simplicity of discussion, we split each qubit on a 
link into two qubits as illustrated in Fig|TJ The string-net 
condensed Z 2 wave function on the original lattice can be 
naturally extended to a state on the new lattice by replac- 
ing a link with 00 and a 1 link with 11. This new state 




FIG. 1: Hexagonal lattice where each link is occupied by two 
qubits. Ideal Z2 model can be defined on this lattice and the 
red, yellow, blue sets of qubits correspond respectively to the 
plaquette, vertex, and link terms in the Z2 Hamiltonian. 



is still an equal weight superposition of all string-nets 
and hence maintains the Z 2 topological order. The new 
system Hamiltonian can be obtained from the old one 
by adding a —ZZ term to each link and expand the pla- 
quette term into a product of X operators on all twelve 
qubits around the plaquette 
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where I denotes all the links and l\ , I2 are the two qubits 
on the link. It is easy to see that the new Hamiltonian in- 
deed has the new string-net condensed state as its ground 
state. The topological order of the system can be de- 
tected from the ground state wave function by calculating 
the topological entanglement entropy of the state. The 
mapping to the new lattice allows this calculation to be 
carried out exactly in a few steps, as illustrated below. 

According to the definition of topological entanglement 
entropy in Ref. l29T or equivalently defined in Ref. l30|) . 
we take out a simply connected region from the whole 
lattice and divide it into three parts A, B, C as shown 
in Fig [l] By calculating the entanglement entropy for 
regions A, B, C, AB, AC, BC, ABC and combining 
them according to 

St P — Sa + Sb + Sc — Sab — Sbc — Sac + Sabc (4) 

we arrive at the topological entanglement entropy S tp of 
the state. The above definition needs to be applied to re- 
gions much larger than the correlation length of the state. 
For the state in consideration, the correlation length is 
zero and the calculation gives the right result for what- 
ever regions we take. We divide the regions by cutting 
through the pair of qubits on boundary links as illus- 
trated in Fig [l] For a region with n outgoing links on 
the boundary, there are 2"~ 1 orthogonal boundary con- 
figurations due to the closed loop constraint of the wave 
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function. Tracing out each boundary configuration con- 
tributes equally and independently to the entropy of the 
region and hence S = n — 1, which includes one term pro- 
portional to the length of the boundary n and one con- 
stant term — 1. The combination in the definition of Stp 
makes sure that the boundary terms of different regions 
cancel out with each other, so topological entanglement 
entropy for the state is then Stp = 1. 




FIG. 2: Tensor product representation of Z2 ground state. 
One tensor is assigned to every three qubits connected to 
the same vertex. Red tensors are on vertices in sublattice 
A and blue ones are in sublattice B. The out-of-plane gray 
links represent the physical indices of the qubits. The tensors 
connect according to the underlying hexagonal lattice. 

This globally entangled state has yet a surprisingly 
simple local representation using the tensor product lan- 
guage. A tensor product state of two dimensional lattice 
model is represented by associating with each lattice site 
m a set of s tensors ITS (a/37...), = 1, 2, s, where s is 
the dimension of local Hilbert space at site m. k is called 
the physical index of the tensor, a/37, the inner indices 
of the tensors, connect to each other and form a graph. 
The wave function (unnormalized) is then given by 
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where C denotes tensor contraction of the inner indices 
according to the connection graph. We omit the inner 
indices here. (We will in most cases ignore normalization 
of wave function in the following discussion and mention 
specifically when normalization is needed.) The tensors 
representing the ground state of the ideal Z2 model can 
be given as follows. We group every three qubits con- 
nected to the same vertex together and assign a rank 
three tensor to each of the eight physical basis states of 
the three qubits. Now every physical index k in Eqn[5] 
is represented with three bits 71177,271.3. The eight tensors 
are: 



T [O oo](000) = l Iron] (Oil) = 1 

T [ ioi ] (101) = l T [110] (110) = l 

all other terms are zero 



(6) 



The physical indices [771772773] correspond to the out-of- 
plane gray links in Fig [2} The inner indices (7172*3) cor- 
respond to the red/blue links in Fig [2] The inner and 
physical indices all have dimension two and are given in 
the same order as shown in Fig[2j Hence the inner indices 
truthfully reflect the configuration of the physical space 
and only configurations with even number of strings at 
each vertex are allowed. It can then be checked that 
only string-net configurations have non-zero amplitude 
in this representation and the amplitude are actually all 
equal. Therefore, the tensors given in [6] indeed represent 
a string-net condensed state-the ideal Z2 ground state. 

This set of tensors serves as a starting point for our 
variational study of topological phase transitions and we 
wish to know what kind of variations of the tensors corre- 
spond to physical perturbations of the Hamiltonian. We 
will study first two specific cases in the following two 
sections. 



B. Z2 model with string tension 

Suppose that we want to know how magnetic field in 
the Z direction might affect topological order. The per- 
turbed Hamiltonian reads: 

H = H Z2 +\J2Zk (7) 

= -£ni-£n^-£^ 2 +A£^ 
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The Zk term commute with the vertex and link term 
rXjeu ^ii ^h^h m t ne unperturbed Hamiltonian, so the 
closed loop constraint is maintained. The ground state 
wave function is still a superposition of string-net config- 
urations, but with different weight. The magnetic field 
adds energy to each string segment, therefore one reason- 
able guess about the ground state is that each string-net 
configuration has weight exponential in its total length 
of string. 
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where the summation is over all string-net configurations 
and L(<pd) is the total string length of a configuration. 
This weighted superposition can still have a simple tensor 
product representation by locally modifying the tensors 
in Eq(6]to 



7[ooo](000) = 9 
T^oi] (101) = 1 
all other terms are zero 



T [ou] (011) = 1 
T [110] (110) = 1 



(9) 



For g > 1, the weight of each string segment is smaller 
by a factor of g^ 1 ^ 2 than that of no string, lowering the 
weight of string-net configurations exponentially. Physi- 
cally, we can imagine this is due to some kind of tension 
in the strings. Therefore, we refer to this wave function 
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as Z2 state with string tension (g). This state cannot be 
the exact ground state of the Hamiltonian given in EqjSJ 
but it is possible that it gives a qualitatively right and 
quantitatively close approximation to the ground state 
and hence might be a good guess for variational study. 
One necessary condition for this conjecture to be true is 
that the topological order of the state remains stable with 
g close to 1, and this is indeed the case as we will show 
below by calculating topological entanglement entropy of 
the state. 
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FIG. 3: Topological entanglement entropy Stp of the Z2 model 
with string tension calculated for the three blue regions as 
shown in Fig [I] 1/g is the weight of each string segment 
relative to vacuum. Stp remains stable for a finite region away 
from the ideal 1.2 TPS and drops sharply to zero at g ~ 1.75. 



In general, this computation is intractable. The equal- 
ity in Eq|4] holds only in the limit of infinitely large re- 
gions A, B, C. Therefore the computation involves di- 
agonalization of exponentially large matrices, each ele- 
ment of which takes exponential time to calculate. For 
Z2 state with string tension, we circumvent this diffi- 
culty by appealing to the special structure of the tensors 
in Eq(9] In order to find the entanglement entropy of 
a region, we map the computation of the eigenvalues of 
the reduced density matrix to the contraction of a two 
dimensional tensor network. While contracting general 
two-dimensional tensor networks is #P hard, the tensor 
networks we are dealing with turn out to be of a spe- 
cial type, called the 'matchgate' tensoi^. 'Matchgate' 
tensors can be contracted efficiently, which leads to an 
efficient algorithm for determining the topological entan- 
glement entropy of this state. The details of the proce- 
dure will be explained in appendix A. 

In the computation, we take the total system size to 
be large enough such that it does not affect the result of 
the computation. Taking the size of regions A, B and 
C to infinity is hard. We manage to carry out the com- 
putation for progressively larger regions as shown in Fig 
[l] with 18, 72 and 162 qubits inside respectively. The 
resulting topological entanglement entropy is plotted in 
Fig[3j We do see a trend of sharper drop of S tp from 1 to 
as the size of the region is increased. At g close to zero 



or g very large when the correlation length is small com- 
pared to the size of the region, the calculated value for 
S tp is reliable and we find that it remains stable within a 
finite range of the ideal Z2 TPS and drops to zero beyond 
certain critical string tension g c . We can decide from the 
plot the critical point g c to be around 1.75. We would 
like to comment that the tensor network representing the 
norm of the Z2 state with string tension g is the same as 
that representing the partition function of classical Ising 
model on triangular lattice with coupling constant J at 
inverse temperature f3 = \ng/{2J). A phase transition 
at g = ^/i is known for this classical model. Our calcu- 
lation for the quantum model confirms this observation 
and shows that the quantum phase transition is in fact 
topological. 

The stability of topological order at g > 1 is a neces- 
sary condition for string tension to correspond to local 
Hamiltonian perturbations. In this particular case, we 
can actually find the corresponding perturbations explic- 
itly. The relation we mentioned above between the Z2 
state with string tension and 2D classical Ising model at 
finite temperature allows the construction of a continu- 
ous family of parent Hamiltonian H{g) for the quantum 
states^!. The Hamiltonian H(g) is local, has the state 
|$f ) as its exact ground state and remains gapped for 
g < g c . Therefore we can conclude that string tension 
can be induced by local perturbations of the Hamilto- 
nian and hence is an allowed variation of the Z2 tensors. 



C. Z2 model with end of strings 

Another simple model one might want to study is the 
Z 2 model with magnetic field perturbation in the X di- 
rection. 

H = H Z2 +\J2x k (10) 

k 

= -£n*-£n^-£ z '^+ A £** 
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The action of the perturbation operator Xk on Z2 ground 
state will flip a link from no string to having a string (or 
back) and hence break one or more closed loops. The per- 
turbed ground state would contain configurations with 
end of strings. In tensor language, this corresponds to 
allowing some odd configurations to be non-zero. Taking 
the translational and rotational symmetries of the Hamil- 
tonian into consideration, one might expect that the fol- 
lowing tensors which assign a small and equal weight e to 
all odd configurations might represent a good trial wave 
function for the ground state. 



7[ooo](000) = 1 
T [101] (101) = l 
T[ooi](001) =6 
T [100] (100) =e 
all others are zero 



T[oii](011) 
7[no](110) 
T[oio](010) 
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Again the inner indices {i^i^j truthfully reflect the con- 
figurations of the physical indices [711712/13]. When e = 0, 
this is reduced to the tensors in the ideal Z 2 TPS. When e 
is non-zero, the wave function contains all possible string 
configurations, closed loop or open string. The weight 
of each string configuration is exponentially small in the 
number of end of strings contained. 



mj = 



(12) 



where the summation is over all possible string configu- 
rations and q{4>sc) is the number of end of strings in a 
particular configuration. 

To see how topological order of the state changes as e 
varies from 0, we again calculate the topological entan- 
glement entropy of the state. In this case, it turned out 
that analytical calculation is possible. The detailed pro- 
cedure is given in appendix B. We find that for any finite 
value of e, when system size goes to infinity, S tp goes to 
zero. Hence topological order is unstable under this kind 
of variation. At first sight this may be a surprising result, 
as we are only changing the tensors locally and we are 
not expected to change the global entanglement pattern 
of the state. However, when we write out the wave func- 
tion explicitly we will see that we have actually induced 
global changes to the state. The wave function in Eq. [12] 
can be expanded in powers of e as 
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where the v's are any vertices in the lattice. |^V 2 '"i) 
is an excited eigenstate of the Z2 Hamiltonian (Eq. [3| 
which minimizes energy of all local terms except the ver- 
tex terms at Vi , Vj and is hence an equal weight superpo- 
sition of all configurations with end of strings at Vi and 
Vj. Note that end of strings always appear in pairs. We 
will call such a pair a defect in the string-net condensate. 
Vi, Vj can be separated by any distance and the number 
of local operations needed to take |$z 2 ) to |<J> Z * 2 ' 3 ) scale 
with this distance. 

On the other hand, with arbitrary local perturbation 
to the dynamics, the Hamiltonian reads 



H — Hj, 



V 



(14) 



where /i u 's are any local operator and 77 is small. The 
perturbed ground state wave function will contain terms 
like I $2*'^) but only with weight ^stance^.v,) _ when 
Vi, Vj are separated by a global distance, the weight will 



be exponentially small. 



Hence a con stan t, finite weight 
is not possible. 
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e for all |$ z y ') as required in Eq. 
Therefore, while we are only modifying the tensors lo- 
cally, we introduce global 'defects' to the state, which 
cannot be the result of any local perturbation to the 
Hamiltonian. We can, of course, design a Hamiltonian 
H e which has \&% 2 ) as its exact ground state using the 
method introduced in Ref. |2U However, H £ will not be 
able to smoothly connect to Hz 2 as e — > 0. 



D. Necessary symmetry condition 

The two kinds of tensor variations we have studied 
have drastically different effects on the topological order 
of the state. While the first type corresponds to local 
perturbations of the Hamiltonian and keeps topological 
order intact, the second type does not have a physical 
correspondence and destroys the topological order com- 
pletely. What leads to such a difference? Given a general 
variation of Z2 tensor, how can we tell if it is allowed? 

We observe that the tensor representing the ideal Z2 
state (Eq[6| has certain inner symmetry, that is, the ten- 
sor is invariant under some non-trivial operations on the 
inner indices, as shown in Fi^4] 




FIG. 4: Symmetry of the Z2 tensor. The tensor representing 
the ideal Z2 state is invariant under the action of Z® Z® Z to 
its inner indices. The variation in string tension (Eq[9| does 
not break this symmetry and topo logi cal order is stable. The 
variation with end of strings (Eq 111 breaks this symmetry 
and destroys topological order 

Z does nothing to the tensor when the index is and 
changes the sign of the tensor when the index is 1. In 
the ideal Z 2 tensor, only even configurations of the in- 
ner indices are non-zero. Hence applying Z at the same 
time to all three inner indices doesn't change the ten- 
sor. That is, Z ® Z ® Z is a symmetry of the tensor. 
As Z <S> Z (g> Z squares to identity, we will say that the 
tensor has Z2 symmetry. Note that we can insert a set 
of unitary operators U, W between any connected links 
in a tensor network without affecting the result of tensor 
contraction and hence the quantity represented by the 
tensor network. Therefore, the Z2 symmetry could take 
any form which is local unitary equivalent to Z ® Z ® Z . 
This Z 2 symmetry is closely related to the closed loop 
constraint of the state. Due to this symmetry, the tensor 
network cannot be 'injective' as defined in Ref. [2H 

Adding string tension to the Z2 tensor (Eq|9| does not 
violate this symmetry, as all the odd terms of inner in- 
dices are still zero. We found that topological order of 
the state is stable with small string tension. On the other 
hand, adding end of strings (Eq |ll[ ) breaks this symme- 
try for any finite e. In general, assume the variation of 
the tensor T contains a Z2 symmetry breaking term dT 
of magnitude S. Such a term would represent an end 
of string in the tensor network. To the leading order in 
8, the wavefunction would contain terms on the order 
0{5 2 ) with dT on two of the sites and T on the others. 
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In the physical space, this would correspond to an open 
string configuration(up to local unitaries at the ends). 
The weight of such a term is 0(5 2 ) even though the two 
sites with dT may be globally apart, hence introducing 
global defects to the wavefunction and breaking topolog- 
ical order. Such defect terms cannot be created by local 
perturbation to the Hamiltonian. Therefore, Z 2 symme- 
try breaking variations to the tensors are not allowed and 
preserving Z 2 symmetry of the tensor is shown to be a 
necessary condition for any variation of the ideal Z 2 ten- 
sor to be physical. This argument is valid for a generic Z 2 
breaking variation. There can be specially designed cases 
where Z 2 breaking variations does not lead to breakdown 
of topological order, e.g. when such variations only occur 
within a finite region of the system or different contri- 
butions to the global defects exactly cancel each other. 
However, for a random Z 2 breaking variation, topologi- 
cal order will be lost and it cannot correspond to local 
perturbation of Hamiltonian. 

The necessity of Z 2 symmetry in the generic case is 
clearly reflected in the following calculation. We ran- 
domly pick tensors in the neighborhood of the ideal Z 2 
tensor and find the topological order of the corresponding 
state numerically. To do this, we make use of a general- 
ization of topological entanglement entropy, the Topolog- 
ical Entanglement Renyi Entropy^. Renyi entropy for a 
reduced density matrix p of order a, where a > 

S a (p) = ^—log[Tr(p a )} (15) 
1 — a 

is a valid measure of entanglement. In the limit of a —> 1, 
it reduces to the usual von Neumann entropy. It was 
shown in Ref. [37] that we can replace von Neumann en- 
tropy with Renyi entropy in the definition of topological 
entanglement entropy (EqQ and still have a valid char- 
acterization of topological order. The resulting quantity, 
topological entanglement Renyi entropy S tpr , does not 
depend on a. We are hence free to choose a for the ease 
of computation and we take it to be 2. The calculation 
of Renyi entropy is mapped to the contraction of a two- 
dimensional tensor network which can be computed ap- 
proximately using the tensor entanglement renormaliza- 
tion algorithm^]. We take the same geometry of regions 
as in Figjl] and the Renyi entropies of different regions 
are then combined in the same way as in Eqj4] to yield 
Stpr- The details of the computation will be described 
in Appendix C. Here we present our result. We restrict 
ourselves to a small neighborhood near the Z 2 tensor 

| r [nin 2 n 3 ](*l*2«3) - T %2 [nin 2 n 3 J (*i*2i 3 )| < 0.1 (16) 

We pick 100 tensors with Z 2 symmetry and plot how 
their topological entanglement Renyi entropy scales with 
reduced region size in the left half of Fig|5] and do the 
same for 100 tensors without Z 2 symmetry in the right 
half of Figj5] We see that for tensors with Z 2 symmetry, 
Stpr approach 1 very quickly as we include more and more 
qubits in the reduced region, while for tensors without Z 2 
symmetry, St pr drops towards as the region gets larger. 



This confirms our statement that Z 2 symmetry is a nec- 
essary condition for any generic variation of Z 2 tensor to 
correspond to physical perturbations of the Hamiltonian 
and hence characterize variations within the topological 
ordered phase. The plot also suggests that Z 2 symmetry 
might be a sufficient condition also. 




reduced region with growing size reduced region with growing size 



FIG. 5: Topological entanglement Renyi entropy (St pr ) calcu- 
lated for three blue regions as shown in Fig[l] Size of regions 
grow from region 1 to 3. The calculation is done for 200 ran- 
dom tensors in the neighborhood of Z2 tensor. 100 of them 
have Z2 symmetry (plotted on the left hand side), while the 
other 100 have not (plotted on the right hand side). For ten- 
sors with Z2 symmetry, Stpr approach 1 very quickly as we 
include more and more qubits in the reduced region, while 
for tensors without Z2 symmetry, Stpr drops towards as the 
region gets larger. 



III. CONCLUSION AND DISCUSSION 

Our result on Z 2 topological order provides useful per- 
spective on the general relation between tensor variation 
and Hamiltonian perturbation. First, it is shown that 
not all variations in tensor correspond to perturbations to 
the Hamiltonian. For the Z 2 model in particular, based 
on our calculation of topological entanglement (Renyi) 
entropy for tensors in the neighborhood of the ideal Z 2 
tensor [see eq. d6b] , we show that one necessary condi- 
tion is that, the tensor is invariant under Z 2 symmetry 
operation Z ® Z ® Z (or any local unitary equivalent op- 
erator) on its inner indices. A generic variation which 
breaks this symmetry cannot be induced by local per- 
turbation of the Hamiltonian and the tensors no longer 
represent state with Z 2 topological order. This gives par- 
tial answer to the question of what kind of variations in 
the Z 2 tensor correspond to physical perturbations to 
the Hamiltonian and hence represent states within the 
same topological ordered phase. Note that we start with 
a particular Hamiltonian in order to better explain the 
property of the state. Our result doesn't depend on this 
particular form of this Hamiltonian and remains valid 
for any local Hamiltonian of the Z 2 topological ordered 
state. (Certain uniformity condition of the Hamiltonian 
must be satisfied, as pointed out in Ref. [2"5"|) Moreover 
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for simplicity of calculation, we restricted ourselves to 
hexagonal lattice in the above discussion. However, the 
Z 2 symmetry requirement is generally true for any lat- 
tice structure and the symmetry operation would take 
the form Z ® Z ® ••• <& Z on all inner indices (or any lo- 
cal unitary equivalent operator). We expect that similar 
necessary symmetry condition also holds for other quan- 
tum double model with gauge symmetry 3 ^. The general- 
ization to other gauge symmetries are discussed in more 
detail in Appendix D. 

This understanding will provide important guidance 
for variational studies of topological order using tensor 
product states. Suppose that, for example, we want to 
find a tensor product state which is the approximate 
ground state of a Hamiltonian with Z 2 topological order. 
It is then very important to search within the set of vari- 
ational tensors that have Z 2 symmetry. If the numerical 
calculation does not carefully preserve this symmetry, we 
might result in a tensor without Z 2 invariance. As the Z 2 
breaking term can be arbitrarily small, the corresponding 
tensor product state might still give good approximation 
to local properties such as energy, but will have totally 
wrong global properties such as topological order. Then 
any attempt to decide the phase diagram based on the 
state would be misleading. 

Finally we would like to note that the symmetry condi- 
tions might not be sufficient. A complete understanding 
of the correspondence between Hamiltonian perturbation 
and tensor variation would be very much desired as it 
might lead to full classification of quantum states and 
quantum phases using the tensor language. 
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Appendix A: Calculating S tp for Z2 model with 
string tension 

In this section, we give detailed procedure of how topo- 
logical entanglement entropy St p of Z 2 model with string 
tension can be calculated using the matchgate tensor 
technique. Following the definition irP3, we take out a re- 
gion (as in Fig [I]) from the hexagonal lattice by breaking 
the to out-going links in half. Due to the closed loop con- 
straint on the wave function, the boundary qubits have 
only 2 m_1 possible configurations a. Regrouping terms 
in the wave function according to different boundary con- 
figurations, we have (up to normalization) 



l$I 2 > = E a *IC 4 M?) 



(17) 



This wave function is automatically in Schmidt- 
decomposition form because for different boundary con- 
figurations |(/>° ut )'s are orthogonal to each other, and 



so are |0™)'s. Knowing the norm and all the a^s would 
enable us to calculate entropy of the reduced density ma- 
trix of the the region. 

Define rank three tensors T, To, Ti with inner dimen- 
sion two as 



T(000) = g 2 
T (000) = g 2 
Tx(lOl) = 1 
all others are 



T(011) = 1 T(101) = 1 T(110) = 1 
To(011) = 1 
Ti(llfJ) = 1 

(18) 

It can be verified that the contraction of T on all vertices 
of the hexagonal lattice gives the norm of | $| ) . To cal- 
culate Q!j for a particular boundary condition Cj, replace 
tensors at the boundary with To if the boundary qubit 
is and with Ti if the qubit is 1 and make sure the first 
inner index is on the boundary link. Contraction of the 
new tensor network will give \tti\ 2 . These three tensors 
satisfy the conditions as defined in Ref. |3S]and are called 
'matchgate' tensors. The contraction of a tensor network 
of N 'matchgate' tensors can be done efficiently (in time 
iV 3 ). Therefore, for a fixed reduced region with boundary 
length to in a system of total size N, the computation 
of entanglement entropy takes time polynomial in N but 
exponential in to. 

We start from a small reduced region (dard blue region 
in Fig |TJ) with a small to, calculate St p and increase the 
total system size N until the change in St p is negligible 
(< 0.01). We repeat this process for different values of g 
and for progressively larger reduced regions (lighter blue 
in Fig[l]). The result is plotted in Fig|S] 



Appendix B: Calculating Stp for 
of strings 



model with end 



Now we show how the calculation of Stp can be carried 
out for Z 2 model with end of strings, analytically. We 
start again with the division of the lattice into sections 
A, B, C as in Fig[T] Without the closed loop constraint, a 
region with to boundary links has 2 m different boundary 
configurations. Rewriting the wave function according to 
different boundary configurations as 



(19) 



we have obtained the Schmidt-decomposed form of the 
wave function and all we need to know to calculate en- 
tropy are the /Vs an d the norm. 

Define rank three tensors S, So, §1 with inner dimen- 
sion two as 



1 



§(000) 
§(001) 

§0(000) = 1 



§(011) 
§(010) 



1 §(101) 
e 2 §(100) 



1 



§(110) 
§(111) 



1 



§o(011) = 1 §o(001) = e 2 §o(010) = e 2 



§i(101) = 1 §i(110) 
all others are 



1 §i(100) =e 2 §i(lll) =e 2 
(20) 



Contraction of tensor S on every vertex of the lattice 
gives the norm of |$| ). To calculate ft for a particu- 
lar boundary condition Cj, replace tensors at the bound- 
ary with So if the boundary qubit is and with Si if 
the qubit is 1 and make sure the first inner index is 
on the boundary link. Contraction of the new tensor 
network will give 1ft | 2 . The contraction of these two- 
dimensional tensor networks can be made efficient by ap- 
plying a Hadamard transformation (|0) — > (|0) + |l))/\/2, 
|1) ->• (|0) - |l))/v / 2)to each of the three inner indices of 
the tensors and transforming them into 



S'(OOO) = \/2(l + e 2 ) §'(111) = 

S o (000) = l + e 2 § (100) 

S (011) = l-e 2 S (lll) 

§i(000) = l + e 2 Si (100) 

Si (Oil) = -l + e 2 Si (111) 
all others are 



72(1- 
-l + e 2 

-1-e 2 

: -1 - 

:l-e 2 



(21) 



It is easy to see that the contraction value of this tensor 
network can be computed analytically, from which we 
know that the entropy of any region with m outgoing 
links is 

l{l + b N *){l + b N °) , {l + b N >){l + b N °) 

s = m -2 iTF * ln " 



1(1-6*0(1- 
2 l + b N 



b N °) . (1 
* In — 



l + b N 
b N <)(l-b N ° 
l + b N 



{22) 



)/(l + e 2 ) and Ni(N ) is the number of 
N = N, + N n is the 



where b = (1 

vertices inside (outside) the region 
total system size. 

Combining the entropy of different regions according 
to Eq. [4] and taking the limit Ni — > oo, N — > oo, we get 
Sl p = whenever e ^ for Z 2 model with end of strings. 



Appendix C: Calculating S tpr for random tensors in 
the neighborhood of Z2 

Redefining topological entanglement entropy in terms 
of Renyi entropy might simplify the calculation. Specif- 
ically, the calculation of Renyi entropy at a — 2 for a 
tensor product state can be mapped to the contraction 
of a single tensor network, which can be computed ef- 
ficiently in one dimension and approximated in two or 
higher dimension. For example, consider a one dimen- 
sional tensor product state (also called a matrix prod- 
uct state) . The lowest dark level in Figj6] gives a side 
view of the state, where red links represent inner indices 
along the one-dimensional chain and gray links repre- 
sent physical indices. Renyi entropy at a = 2 is de- 
fined as S 2 (p) = -log[Tr{p 2 )}. To find out Tr(p 2 ), 
we stack four copies of the states together as in FigjBJ 
connect corresponding physical indices outside the re- 
duced region between levels 1&2, 3&4 and connect those 
within the reduced region between levels 1&4, 2&3. Con- 
traction of this four-layer one-dimensional tensor net- 
work gives Tr{p 2 ). For two dimensional tensor product 




FIG. 6: Tensor network for calculating Renyi entropy (a = 2) 
of the big site in a one-dimensional tensor product state. 
The contraction of the tensor network gives Tr(p 2 ) = 
exp(—S-2(p)), where p is the reduced density matrix . The 
lowest level represents the state, where red links represent 
inner indices along the one-dimensional chain and gray links 
represent physical indices. Four copies of the state are stacked 
together and corresponding physical indices are connected be- 
tween the levels. For physical indices outside the reduced re- 
gion, the connection is between levels 1&2 and levels 3&4. For 
those in the reduced region, the connection is between levels 
1&4 and levels 2&3. 



states, the generalization is straightforward. The only 
difference is that now we have to contract a four-layer 
two-dimensional tensor network. To this end, we apply 
the Tensor Entanglement Renormalization Algorithm^. 
Having obtained the Renyi entropy for different regions, 
we then combine them to get S tpr . 



Appendix D: Gauge symmetry of tensor product 
states and topological order 

This section discusses in general the relation between 
gauge symmetries of tensor product states and topologi- 
cal order, and the implication of our result on other topo- 
logical ordered models with gauge symmetry. 

For a tensor product state, the network of tensors 
which represents the same state is not unique. In partic- 
ular, if we change a pair of connected tensors by rotating 
the basis of one of the connected inner index with an 
invertible operator A and rotating the other connected 
inner index with operator A -1 , any tensor trace would 
remain unchanged and hence the tensor product state re- 
mains the same. This corresponds to inserting a pair of 
invertible operators A, A" 1 onto any link in the graphi- 
cal representation of the state. Following the definition in 
Ref. |3H1 this is called a gauge transformations of the ten- 
sor product state, which form a very large group. Hence 
the correspondence between the tensor network and the 
physical state is many-to-one. As a result, the variation 
energy as a function of tensors has a very large symme- 
try: the variation energy is invariant under the gauge 
transformations. 

On the other hand, when we try to find the best de- 
scription of ground state for a model Hamiltonian by min- 
imizing energy with respect to the variations of the ten- 
sors, the tensors that minimize the average energy may 
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not be invariant under all the gauge transformations and 
in general have much less symmetry. For example, in the 
ideal Z 2 case, the tensors are only invariant if we insert Z, 
Z~ x to all the links in the 2-dimensional graph. General- 
izing this to any symmetry group and to any dimension 
d, we define the d-dimensional Invariant Gauge Group (d- 
IGG). The d-IGG is nothing but the invariant group of 
the tensors under gauge transformations. Thus the min- 
imization of the average energy leads to a spontaneous 
symmetry breaking. The d-IGG's are the unbroken sym- 
metry of the tensors that describe the ground state. As 
we change the Hamiltonian, the tensors that minimize 
the average energy may have some different symmetry 
structures described by different d-IGG's. As is shown 
in Ref. 1381 when d equals the dimension of system space 
dspace, d space -lGG (such as the Z2 symmetry discussed 
in this paper) can be used to determine the topologi- 
cal orders of a tensor product state. A closely related 
concept is discussed in Ref. Therefore, a change in 
d space -lGG will in general represent a change in topo- 
logical order. Apart from d space TGG, the tensors might 
also have lower dimensional IGG's. For example, if we 
trivially map every inner index i to ii(i — 0, 1) in the Z2 
tensor, the tensors still represent the same state but have 
a 0-IGG ZZ in additional to its 2-IGG. However, we be- 
lieve that such lower dimensional IGG's are not related 
to the topological order in two dimension and breaking 
them will not lead to a change in topological order. 



Note that in order to use <i sapce TGG of a tensor net- 
work to decide topological order, we only require that 
the network is composed of patches of tensors which are 
invariant under certain gauge transformations. It is not 
necessary that every single tensor is d sapce -lGG invari- 
ant. However, in the generic case, if the single tensors 
do not have special symmetry structure, it is not possi- 
ble to have d space -IGG invariance on a bigger patch. As 
discussed in Ref. [2H such a tensor network will gener- 
ically satisfy a condition called 'injectivity', i.e. for a 
large enough region in the network, when the single ten- 
sors are contracted together to form a new tensor, the 
set of tensor vectors labeled by their physical indices will 
span the full tensor space of the n outgoing inner indices 
of the region. Therefore, the tensor network cannot have 
nontrivial c? spQce -IGG. In order for a bigger patch in the 
network to have d space TGG invariance, it is in general 
necessary for every tensor to be d space TGG invariant. 

Hence, we believe that the invariance of every tensor 
under <i space -IGG is a more general necessary conditions 
for generic variations of the tensor to correspond to phys- 
ical perturbations of the Hamiltonian. Breaking of the 
dspace-IGG invariance of the tensors will in general cor- 
respond to a change in topological order. Therefore in a 
numerical variational calculation it is very important to 
preserve the d space -IGG invariance. Otherwise we would 
not be able to correctly determine the topological order 
of the resulting state from the tensors. 
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